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Abstract 

We investigate dynamical symmetry breaking of the Gross-Neveu model in the light- 
front formalism without introducing auxiliary fields. While this system cannot have zero- 
mode constraints, we find that a nontrivial solution to the constraint on nondynamical 
spinor fields is responsible for symmetry breaking. The fermionic constraint is solved by 
systematic 1/N expansion using the boson expansion method as a technique. Carefully 
treating the infrared divergence, we obtain a nonzero vacuum expectation value for fermion 
condensate in the leading order. We derive the 't Hooft equation including the effect of 
condensation, and determine the Hamiltonian consistently with the equation of motion. 
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1 Introduction 



One of the most crucial problems in the light-front (LF) formalism is how to describe spon- 
taneously symmetry breaking on the trivial vacuum. The idea that solving a constraint 
on the longitudinal zero mode (zero-mode constraint) may provide us with mechanism for 
it [1|] has been examined by several authors in 1+1 dimensional scalar models 0, [|. If the 
zero mode has a nontrivial c-number part, it gives a nonzero vacuum expectation value of 
fields. Although the nonperturbative calculation is still difficult in 3+1 dimensions, this 
approach is indispensable for comparison with the other approach [5| where , instead 
of solving zero mode constraints, nontrivial vacuum effects are supposed to be able to 
be incorporated by suitable counter terms with symmetry consideration. Then, from the 
standpoint of the zero-mode approach, a straightforward question arises: if there exists 
no bosonic field in the system, how one can describe the symmetry breaking even without 
zero-mode constraints. This question motivated us to investigate the simplest fermionic 
model, the (massive) Gross-Neveu model with global U(iV) symmetry ||: 



,2 



c = y(iYd„-m )v + ^(my. (i.i) 

In the equal-time formulation, we know that when the bare mass is absent, the discrete 
symmetry \1/ — > 7 5 \I/ breaks dynamically via (tyty) ^ 0. So far the study of the light- 
front Gross-Neveu model has been a little controversial. First, the same phenomenon 
was successfully observed using the light-like quantization surface which approaches the 
light-cone surface as a limit |7], EJ. On the other hand, exactly on the light-front, it has 
been an unsettled problem whether we can obtain a nonzero condensate or not. as well 



as the renormalization procedure || 10]. However, Maedan [11, 12] obtained the same 
result as that of the equal-time formulation. He solved the zero mode constraint, which 
was enabled by the introduction of a bosonic auxiliary field. 

In this letter, we shall again make an analysis of the Gross-Neveu model without 
auxiliary fields included. To compare with the previous works and to see carefully the 
way condensation emerges, we put the system in a box of length L and impose antiperiodic 
boundary condition. 
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2 Formulation of the problem 

The Gross-Neveu model on the LF has two features characteristic of the four-fermi inter- 
action. In the LF formalism, half of the spinor field is a dependent variable and should be 
constrained. The constraint in the Gross-Neveu model, however, is difficult to solve due 
to the nonlinearity of the four-fermi interaction. Indeed, the Euler-Lagrange equations 
are given by 

id+il> = + (2.1) 

id ~ x = te-y^x + xty)}^, (2.2) 

where \l/ a = 2 _1//4 (^ a , Xa) T , (a = 1, • • • , N) and the light-front derivatives are d± = 
d/dx ± ,x ± = (x°±a; 1 )/v / 2, x~ G [—L,L], and the 7 matrices 70 = 01,71 = £03. Equation 
( p.2|) includes only spatial derivative and thus is a constraint, which we call hereafter as 
a "fermionic constraint". 

The Gross-Neveu model can also be defined as a limit of the Yukawa-like theory |12| . 



At the limit, the scalar field of the system becomes an auxiliary field. There exist two 
constraints. One is the zero-mode constraint for the auxiliary field and the other is the 
fermionic constraint which is linear in terms of the spinor fields. While the fermionic 
constraint is easily solved, the zero-mode constraint is complicated and gives a nonzero 
value of condensation. However, the solution to the fermionic constraint involves the zero- 
mode of the auxiliary field and, thus in a strict sense, it is not solved at the first step. 
What was done in |CJ was to eliminate the nondynamical degrees of freedom from the 
coupled equations by a two-step approach. Therefore we expect that, without auxiliary 
fields, the fermionic constraint must be responsible for dynamical symmetry breaking and 
so we shall carefully solve eq. (|2.2j). 

Another problem is vanishing of the classical Hamiltonian in the massless limit. Ac- 
cording to Dirac's procedure, the Hamiltonian on the constrained surface is given by 

H = ^= fdx-Wx + xty), (2-3) 



2V2 



L 



where x should be substituted by a solution of ( |2.2|) . Thus, the Hamiltonian vanishes 



in the massless limit |10[. However, if the spontaneous symmetry breaking occurs, the 



fermion acquires a nonzero mass and even if the bare mass goes to zero, the massless 
limit must have nontrivial Hamiltonian. Nevertheless, it would be difficult to observe 
such effects in the classical Hamiltonian since the symmetry breaking in the Gross-Neveu 
model takes place only in the quantum level. Therefore we decide to start from the Euler- 
Lagrange equations as in Ref . || . The Hamiltonian will be constructed consistently with 
the equation of motion after we solve the fermionic constraint. 

Since the fermion condensate is given as the vacuum expectation value of the fermion 
bilinear operator, it would be better to rewrite ( |2.1| ) and ( |2.2| ) in terms of bifermion 
operators. We introduce U(iV) singlet bilocal operators at the equal light-front time: 

N 

M(x-,y-;x + ) = £ ^(aT, x+)^ a {y~, x + ), (2.4) 

0=1 

N 

C(x-,y-;x + ) = ^2^ a (x-,x + ) Xa (y-,x + ). (2.5) 

a=l 

We also introduce T(x,y) = C(x,y) + C*(x,y) for convenience. From now on, we work 
at the fixed light-front time and omit x + -dependence and superscripts for the light-front 
spatial variables unless needed. 

The equations for the bilocal operators are written as 

2 

id + M(x, y) = ^= (C{x, y) - C% x)) - 9 - (C(x, y)T(y, y) - T(x, x)C\y, x)) , (2.6) 

i-^T(x,y) = ^L(M(x,y)-M(y,x)) - 9 ^(M{x,y)T{y,y) -T(y,y)M(y,x)). (2.7) 

We define the theory by these equations with this ordering and a quantization condition 
on the dynamical fermion: 

{ip a (x), ii{y)}x+= y + = Sab5(x~ - y~), {^pa(x) , ip b (y)} x + =y + = 0. (2.8) 

Fourier expansions of the bifermion operators are defined by 

M(x,y) = — e- iknX e- ikmy M(n,m), (2.9) 

~ i 

and so on, where k+ = irn/L and n, m are half integers due to the anti-periodic boundary 
condition ip(L) = —ip(—L). Momentum representation of the above equations are 

id + M(n,m) = ^-^(C(n,m) — C\—m, —nfj 

2 

- Jf E { C ( n > m-k- l)T{k, I) - T(k, /)C f (-m, k + I - n)) ,(2.10) 

k,i 
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and 

k m T(n,m) = ^-^(^M(n,m) — M(m,n)j 



m 
AL 



(M{n,m-k-l)T{k,l)-T{k,l)M{m-k-l,n)){2.11) 
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3 Solution to the fermion constraint 

We solve the fermionic constraint ( p.llj ) by using the 1/N expansion. Let us expand the 
bilocal operators as 

M(n,m) = N^[-=] fi (p \n,m), (3.1) 



p =o \ v N , 



and similarly, C(n, m) and T{n, m) are expanded by c^(n, m) and t^(n, m), respectively. 
Expansion of M(n, m) can be given by the boson expansion method which is a familiar 
technique in the many-body physics [13]. Among various ways of the boson expansions, 
the Holstein-Primakoff type is the most useful for large N theories. Introducing bosonic 
operators 

B{m,n 2 ), B\m x ,m 2 ) = 6 niimi 6 n2 , m2 , B{m,n 2 ), B{mx,m 2 ) =0, (3.2) 
M(n, m)'s are represented as follows: 

: M_ + ( ni ,n 2 ) : = £ B\-m, k)B(n 2 , k) = A(n 2 , -m), (3.3) 

1,-1 2 ... 
K — 2'2' 

:M + _(m,n 2 ): = - £ B\k,-n 2 )B{k,n x ), (3.4) 

K- 2>2> 

:M ++ (m,n 2 ): = £ (v / iV^4)(^2, k)B(k,n 1 ), (3.5) 

K — 2'2' 

:M__Kn 2 ): = 2 5 t (^-«2)(v / iV 3 TA)(A;,-n 1 ), (3.6) 



fc ~2'2- 



where the suffices imply the sign of the momentum and the normal order is defined on the 
Fock vacuum. The right-hand-sides of these are determined so that this representation 
satisfies the algebra of : M(n,m) :. If we expand eqs. (|3.5|) and ( |3.6|) in terms of 1/N, 
each order is explicitly given by B and B' unlike the naive 1/N expansion of M(n, m) in 
Ref. [[J. Therefore, we can express also in terms of the bosonic operators B and B^ 
in principle. 



Using //(°)(n, m) = 9(n)5 n+m)0 , the lowest order of the fermion constraint is given by 
k m t^\n,m) = ^e(n)6 n+mfi -^re(n) £ t®(k,n + m-k), (3.7) 

where #g = g 2 N and e(n) = 0(n) — 9(—n). We find that t^(n, m) is a c-number because 
there are no operators in this equation. The c-number part of a bilocal operator is a 
function of only x — y due to the translational invariance of the vacuum, which implies 
t(°)(n,m) =t(°\n, -n)6 n+mfi . 

By the way, if the lowest order equation has a nontrivial solution, it gives rise to the 
physical fermion mass 

M = m -g 2 (m)=m --^- £ t®(k,-k). (3.8) 

2 ^ 2L kez+i 

Inserting this and m = —n into (|3.7| ), and summing over n G Z + 1/2, we obtain an 
equation 

1 1 M „ Ak _ 

~gJ~2^M-m _fs (3 ' 9) 

2'2'"" 

where A A; = n/L. If we set m = naively, this equation seems to become independent 
of M and thus we cannot determine M by this equation. Or if we solve ( |3.9|) in terms 
of M, and take the m — * limit, M also seems to vanish. However, this observation is 
not correct because this equation is not well-defined until the summation is regularized. 
Indeed, by carefully treating the divergences, this equation gives a nonzero value for M 
even in the mo = case. 

Since eq. (|3.9| ) is independent of the box length L, we can evaluate it in any box 



size. Nevertheless we evaluate it in the continuum limit L — > oo because it is what we 
want eventually. Then, we find that the summation has both infrared and ultraviolet 
divergences. We regularize it as follows: 

Ak _ 1 rA dk OA 2 

L r = L ^ M /o — > L — = ln TH- ( 3 - 10 ) 

n- 



,,s_...K ^ m + l/2 JmI k M 2 



2 ' 2 ' 



This is considered to be the parity invariant regular ization. Since the parity trans- 
formation is the exchange of k- and k + , the cut off < A should be paired with 
k + = M 2 /2/c_ < A, which inevitably relates the UV and IR cutoffs. Here the dispersion 
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of the fermion is thought to be 2k + k_ = M 2 because the fermion acquires the physical 
mass ( |3.8| ). This procedure corresponds to imposing a consistency condition. If we in- 
troduce UV and IR cutoffs independently, it generally breaks the parity invariance and 
we cannot yield the correct result. It should be commented that eq. (|3.9|) in the mass- 
less case is identical with the leading order of the zero mode constraint in the auxiliary 
field approach [Q. There, the infrared regularization was performed by the heat kernel 



method, which gives the same result. 

Renormalization is implemented as follows. The coupling constant is renormalized as 

1 1 1 , 2A 2 1 . 

SM = ^-^V + -' (3 ' u) 

/ 1 2A 2 l\ 

9t = Z{n)gl{p), Z -1(0) = ! + ^ /_]*—- -J, (3.12) 

which shows the asymptotic freedom. Mass can be renormalized by the same factor: 
m = Z(/i)m^(/i). Rewriting eq. ( |3.9|) in terms of the renormalized quantities, we obtain 
an equation for the physical fermion mass; 

M 2 _ M — m R 2n 

In the rriR = case, the solution is 



l ^ = 2 -^r JL ^- ( 3 - 13 ) 



M = /iexp J 1 - 4" \ = M . (3.14) 
I 9r) 

This is the same result as that of the original work by Gross and Neveu |§. In the massive 



case, using the renormalization invariant parameter M , eq. ( 3.13 ) can be written as 

Mf, g\M 

This equation has a nonvanishing solution. If m R /M <§C 1, it is M ~ M + (n / ' g\)mn 
which smoothly goes to M as m R — > 0. Eventually, is determined from eq. (|3.7|) as 



m , . M 6{n) r 

V2 fc n 



Similarly, the fermionic constraint in the next leading order 
k m t^ (n, ra) = — -=(^^(71, m) — //^(m, n)) 



V2' 

^-e(n) t^(k,n + m-k), (3.17) 
feez+l 



is solved as 

c^(n, m) = —=.- — n^ x \n, m) (3.1 

V 2 k m 

where K = n + m. This can be expressed in terms of B(n, m) and B>(n, m) by using 

m) = B(m, n), ^„(n, m) — B\— n, — m), (3.19) 

^l(n,m) = 0, /z^(n,m) = 0. (3.20) 

Note that there is no mass correction in this order due to fr- 1 ' (n, — n) = 0. Also we 
can easily obtain the higher orders because the equation is always linear in the highest 
order c^ p \ Thus C(n,m) can be represented only by the bosonic variables B(n,m) and 
B^(n, m). 

4 Hamiltonian and the 't Hooft Equation 

We explicitly solved the fermionic constraint up to the next leading order and C(n,m) 

can be represented by B and B\ The next work is to rewrite the equation of motion in 

terms of the dynamical variables and construct the Hamiltonian from it. 

The nontrivial leading contribution of the equation of motion ( |2.10|) is 

• n m/ \ M 2 f I 1\ m . . 
id+p'fam) = — I — + — J fi { '{n, m) 

A 

4L y/2 vv "' " y \k n k 

where 

\ 0<n<\K\ n 

x^= J: T^U l \n,K-n)-^\K-n,n)). (4.2) 
V2 nez+ i k K . n 

Rewriting this in terms of the bosonic variables, we obtain equations of motion for B(n, m) 
and B\n, m): 



zd + B(n,m) = J2 K kfB(k,l), (4.3) 

k,l,n,m>0 

id + B\n,m) = - E KpB\k,l), (4.4) 



k,l,n,m>0 



where the matrix is 



K n,m M 2 / 1 1 

K k,i = ~y [tt + T~ i 5n ' k6ri 



m R gl AA;\%|M 2 /1 lWl 1\ 

^T + ^T 7— TT - ^ - I - / \ T. 17.) °k+i,n+m- (4.5) 



M 2lT Q< ^ Kl k m J 4L 2 V&„ k m J \k k ki 
It is easy to find the Hamiltonian which gives these as the Heisenberg equations of motions: 

H= Yl 5 t (n,m)K"; / m 5(A;,/) + 0(Ar- 1 / 2 ). (4.6) 

n,m,k,l>0 

Let us introduce the collective operators 

&£>t = £ & a \n)B\n,K-n), (a = l,...,K), (4.7) 

o<n<i<: 

where the wave function 3>w ( ra ) i s assumed to satisfy the orthogonality and the complete- 
ness. Then bx satisfy the bosonic relations [b^\ b^,^} = S aja '8K,K'- Note that bx does 
not have the zero mode because both of the arguments of B should be positive. The wave 
functions are determined so as to diagonalize the Hamiltonian. The eigenvalue equation 
gives the 't Hooft equation. If we introduce rescaled variables x n = k n /kK for < n < K, 
we obtain 

2P + k K $(n) = M 2 ( — + — - — ) $(n) (4. 



2tx Q <m<K x m J Air \x n 1 — x n J 0< i <K \xi 1 — xiJ 

where we take the time dependence of B^ to be e lP+x+ . Taking the continuum limit 
lim^i^oo nK/L = P_, the 't Hooft equation becomes 

2P + P_0(x) = M 2 (— + ——] <j>(x) (4.9) 

\x 1 — X ) 

t'^R ii p^y 1 j|m^ / 1 — 1_\ ri / 1 — 1_\ 

\M 2ix h x) 4vr \x 1-xJJo y \y l-y) Vyy! 

One can easily observe that when x = 1/2, the second term of the 't Hooft equation 
vanish. Thus the invariant mass of this state is 2M which is known to exist as Gross and 
Neveu discussed ||. The wave function of this state is <f>(x) = S(x — \) and far from the 
collective one. In this order, the state 6^|0> = B\K/2,K/2)\0) can be easily translated 
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into a bifermion state by using eq. ( 3.19Q , which means that the state can be understood 
as the constituent state. 

The result ( |4.9|) is different from that of || by the divergent factor of the second term, 
which, however, does not affect the above solution because of the vanishing of the second 
term. This discrepancy originates from the different renormalization procedures. If we 
renormalize only the UV divergence, we do not have such a divergent factor ||. The 
divergent factor seems to be inevitable in our renormalization prescription. However we 
expect that it cancels with the possible infrared divergence of the last integral. Thus 
our treatment is different from the usual way of the 't Hooft equation. Our expectation 
may be strengthened by a little more nontrivial example in the chiral Gross-Neveu model 
Ant = 5 ,2 /2[(^ ,x I / ) 2 — (\I/7 5 \I/) 2 ]. There is the chiral symmetry when uiq = and if it breaks 
spontaneously, there appears a Nambu-Goldstone boson. Although there is no NG boson 



in two-dimensional field theories |fl4| , y|, we can observe it in the leading order. In the 



same way as that of the Gross-Neveu model, the 't Hooft equation becomes 

2P + P_0(x) = M 2 (— + - ) <f>{x) (4.10) 



x 1 — X 



'm R g% r 1 dx\ 1 g\M 2 r 1 ( 1 1 
~M + 2Wo ~) 2tt Jo V [ xu + (1 - x)(l 



<t>{v)- 



\xy (l-x){l-y) / 

Although there emerges in the second term the same factor as in eq. (|4.9f) , we find 
that 4>(x) = constant is the exact solution for massless case. The invariant mass of this 
state is zero and this solution corresponds to the Nambu-Goldstone boson. Precisely the 
divergence of the last integral cancels with the divergent factor. 

5 Discussion 

We calculated the nonzero value of condensation in the Gross-Neveu model on the light 
front without introducing auxiliary fields. We found that the nontrivial solution to the 
fermionic constraint lead to the fermion condensation. In solving the fermionic constraint 
we used the boson expansion method for the systematic 1/N expansion, although it was 
not essential for the analysis up to the next leading order. 

To obtain the nonzero value of condensation, it is very important to treat the infrared 
divergence with great care. It seems suggestive that the prescription for the infrared 
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divergence can be set to be related with symmetry such as the parity invariance. Such 
symmetry consideration will be indispensable to obtain the correct value of the conden- 
sation, which will be true for another approach without zero-modes [f|, ||. 

The emergence of the nonzero condensate is connected with the renormalization pre- 
scription and affects on the 't Hooft equation. That is, the equation cannot help hav- 
ing a divergent factor. However, the 't Hooft equations we obtained has not only a 
well known solution in the Gross-Neveu model but also a solution corresponding to the 
Nambu-Goldstone mode in the chiral Gross-Neveu model. Usually, the 't Hooft equation 
is derived without effects of condensation included and the fermion condensate is obtained 



by, for example, the sum- rule calculation [16|]. It would be interesting to investigate the 
relation between two approaches. 

We constructed the Hamiltonian consistently with the equation of motion after solving 
the fermionic constraints. Therefor we cannot calculate the vacuum energy and cannot 
determine which vacuum is realized. Indeed, there is also a symmetric solution, M = 
in the rriR = limit. However, in this case, if we start from the very massive fermion and 
decrease the mass, we can uniquely determine the vacuum and it is the broken phase. 
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